Timelike 7*^ — > A form factors and A Dalitz decay 



o 

(N 



G. Ramalho 1 and M. T. Pena 1 ' 2 

1 CFTP, Instituto Superior Tecnico, Universidade Tecnica de Lisboa, 
Av. Rovisco Pais, 1049-001 Lisboa, Portugal and 
2 Physics Department, Instituto Superior Tecnico, 
Universidade Tecnica de Lisboa, Av. Rovisco Pais, 1049-001 Lisboa, Portugal 

(Dated: March 3, 2013) 

We extend a covariant model, tested before in the spacelike region for the physical and lattice 
QCD regimes, to a calculation of the 7* AT — > A reaction in the timelike region, where the square of 
the transfered momentum, q 2 , is positive (q 2 > 0). We estimate the Dalitz decay A — > Ne + e~ and 
the A distribution mass distribution function. The results presented here can be used to simulate 
the NN —¥ NNe + e~ reactions at moderate beam kinetic energies. 
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I. INTRODUCTION 

Electromagnetic reactions which induce excited states 
of the nucleon are important tools to study hadron struc- 
ture, and define an intense activity at modern accelera- 
tor facilities, namely at MAMI, MIT/Bates and Jeffer- 
son Lab. An enormous progress in these experimental 
studies has been achieved in the recent years, leading to 
very accurate sets of data on 7*iV — > N* excitation re- 
actions, for several N* resonances at low and high Q 2 
(with Q 2 = -q 2 ) [lHl. This wealth of new experi- 
mental data establishes new challenges for the theoret- 
ical models and calculations, since, in the impossibility 
of solving exact QCD in the momentum transfer regime 
of Q 2 = — 10 GeV 2 , reliable effective and phenomeno- 
logical approaches are unavoidable. 

In this context, we developed a covariant constituent 
quark model for the baryons within the spectator frame- 
work Q for a quark-diquark system . Because the 
construction of the electromagnetic current is based upon 
the vector meson dominance mechanism, it was possible 
to apply the model also to the lattice QCD regime in a 
domain of unphysical large pion masses [1, [HI, Tl3| . The 
model is constrained by the physical data for the nucleon 
and 7*7V -> A datafp, 0, 03 as well as the j*N ^ A 
lattice QCD data |12fl . The evidence of the predictive 
power of the model comes from its results obtained with- 
out further parameter tunin g, fo r the form factors of the 
reactions 7* A — > A(1600) 14| as well as the reaction 
7* TV — > N* where N* can be first radial excitation of 
the nucleon AT*(1440), and the negative parity partner 
of the nucleon AT* (1535) [H-[l3- Moreover, the exten- 
sion of the model to the strangeness sector was successful 
in the description of the baryon octet [lH and baryon de- 
cuplet form factors [8| . All parameters in the model have 
a straightforward interpretation: they give, for instance, 
the momentum scales that determine the extension of 
the particle, and the coupling of the photon with the 
constituent quark. 

Importantly, the information extracted from the elec- 



tromagnetic excitation reactions is also relevant for the 
interpretation of production processes induced by strong 
probes. Of particular interest is the study of NN colli- 
sions in elementary nucleon- nucleon reactions and in the 
nuclear medium [19- 26]. In this sector, the HADES ex- 
periments of heavy-ion collisions in the 1-2 GeV range 
play a unique role in accessing nuclear medium modifi- 
cations at intermediate and high energies [26T - [28j . Fur- 
thermore, in the near future, the FAIR facility will 
expand these experiments further to a higher energy 
regime [lij [U [23, [24| . In both cases, independently of 
the energy domain under scrutiny, the di-lcpton chan- 
nel, the first of which is the low mass di-electron chan- 
nel, is one of the interesting production channels from 
heavy-ion collisions expected to signal in-medium behav- 
ior. It is crucial, for the interpretation of both present 
and planned di-lepton production data from heavy-ion 
experiments at intermediate energies, to have a reliable 
baseline made of experimental reference from nucleon- 
nucleon scattering - one of the objectives of the HADES 
experiments. 

Nonetheless, one also needs an extension of the knowl- 
edge gained from the experimental studies on elementary 
electromagnetic transitions, to the timelike region (Q 2 = 
—q 2 < 0) since the description of the NN — > NNe + e~ 
reaction involves baryon electromagnetic transition form 
factors in that kinematic regime [2l|, [24[ • Natural re- 
quirements of this extension is that it is well constrained, 
appears to be robust when tested by its predictions, and 
allows a direct physical interpretation of the parameters 
involved. 

Therefore, in this work we extend our form factor cal- 
culations for the 7* A — > A reaction to the timclikc re- 
gion and calculate the partial width decays A — > -fN and 
A — » e + e~ N. We follow the procedure of the standard 
simulation packages that treat the low mass di-electron 
production data as a Dalitz decay following a resonance 
excitation [ijj]. We note, in particular, that according to 
Ref. [l!| the fraction of di-lepton events compared to the 
hadronic channels depends significantly on the resonance 
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mass W, and on the details of the dependence on q 2 of the 
transition form factors, two features that call for studies 
as the one we describe here, where such sensitivities are 
investigated. 

We start with the valence quark model presented in 
Ref. [IH for the j*N — ¥ A reaction. That model has 
two important ingredients: the contributions from the 
quark core, and a contribution of the pion cloud dress- 
ing. In the quark core component, the A system has 
a quark-diquark effective structure with an S-wave or- 
bital state and small D-wave admixtures [H, 0, E3] • We 
take here only the dominant S-state contribution which 
is largely responsible for the magnetic dipole transition 
form factor G* M , because D-wave states give only small 
contributions to the A wave function (< 1%) Since 
the electric and Coulomb quadrupole form factors are de- 
termined by the D-state admixture coefficients [3] , in the 
S-state approximation those two sub-leading form factors 
become identically zero. This is a reasonable approxima- 
tion because they are indeed small when compared with 
G*m ■ To the quark core contributions it is neces- 
sary to add contributions from the pion cloud, in order 
to describe the reaction in the physical regime for small 
momentum transfer 0, @] ■ An important feature of our 
model is that it describes the 7*iV — > A reaction in the 
physical and lattice QCD regimes. 

The extension of the model to the timelike region pre- 
sented here is done directly by extrapolating the valence 
quark model @ , fixed in the spacelike region, to the kine- 
matic conditions of the timelike region. This means that 
an arbitrary mass W of the A is taken to replace its phys- 
ical mass value. We also need to generalize the photon- 
quark current to the timelike region, while keeping its 
vector meson dominance parametrization [f| Q . This is 
done by adding a finite width to the vector meson pole 
of the current. As for the pion cloud contributions, we 
study two different extensions to the timelike region. Al- 
though similar in the spacelike regime @, 0] , they have 
very different behaviors in the timelike region. From the 
obtained results we conclude that the model which in- 
cludes the %PT constraints is favored. 

This work is organized in the following way: in Sec. |TT] 
we introduce the formalism that relates the A Dalitz de- 
cay with the -f*N — > A electromagnetic form factors; in 
Sec. IHII the spectator quark model is introduced and the 
explicit expressions for the form factors are presented; 
in Sec. IIVI we show our results for the decay widths of 
A — » 7./V and A — > e + e~N, and for the A mass distribu- 
tion, as a function of W; finally in Sec.|V]we summarize 
and draw our conclusions. 



II. BREIT-WIGNER DISTRIBUTION FOR THE 

A RESONANCE 

In the simulations of NN reactions one has to take into 
account the intermediate excitations of the nucleon, and 
the A (spin and isospin 3/2) resonance is the first relevant 



one [19H221 124 l26j . For that purpose we calculated the 
contribution of the A(1232) state to the cross section, for 
an arbitrary resonance mass W which can differ from the 
resonance pole (defining the mass Ma). The most usual 
ansatz is the relativistic Breit-Wigner distribution [l9], 
USUI, given by 



9a(W)=A 



W 2 Y tot (W) 



(W 2 - Mi) 2 + W 2 [T tot (W)Y 



(2.1) 



where T to t is the total width, dependent of W, and A 
is a normalization factor determined by the condition 
/ dWg&(W) — 1. The total width can be decomposed 
into the contributions from the independent decay chan- 
nels pjj: 

T tot (W) = T nN (W) + T jN {W) + T e+e - N (W), (2.2) 

respectively for the decays A — > irN, A — > 7 TV (7 repre- 
sents a real photon) and A — > e + e~N. 

The dominant process is the decay A — > ttN, which 
can be described by the well known ansatz [l9|, [3l[ 



M A ( q*(W) V ( u(W) x 2 



W \qn(M A )J \v(M A ) 



(2.3) 



where q n (W) is the pion momentum for the decay of a A 
with mass W, and T® N is the A — > ttN partial width for 
the physical A [r^ = T wN (M A )]- The function v(W) 
is a phenomenological function given by 



v(W) 



rf 2 



(2.4) 



where j3 as a cutoff parameter. Following Refs. [19|, |31( 
we use $ = 300 MeV. 

As for the components T 1 n{W) and F e + e - N (W), they 
will be determined by the A Dalitz decay, as described 
next. 



A. A Dalitz decay 

The A Dalitz decay can be expressed in terms of the 
function r 7 *jv(g; W), where 7*jV is a short notation for 
the reaction A — ¥ 7*iV, and 7* represents a virtual pho- 
ton, with squared momentum q 2 > (i.e. timelike). The 
variable q is defined by q = yfq 2 ■ The case q 2 — corre- 
sponds to the real photon limit. 

The r 7 *jv(g; W) function can be written QjlHH as 

Y,, N { Tl W) = ^ W ^ 2 ^y^y^Gr^-W)] 2 , 



where M is the nucleon mass, a 
structure constant, and 



y± = (W±M) 2 



(2.5) 

1/137 the fine- 



(2.6) 
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The function \Gr(q 2 ] W)\ depends on the 7* A -> A tran- 
sition form factors: G* M (magnetic dipole), G%> (electric 
quadrupole) and G* c (Coulomb quadrupole) [33|, and is 
given by 

\G T (q 2 ;W)\ 2 = \G* M (q 2 ;W)\ 2 +3\G* E (q 2 ;W)\ 2 

„2 



2W 2 



G* c (q z ;W) 



(2.7) 



In this equation we note that the contribution of each 
form factor will always be real and positive, even if the 
form factors are complex. 

Equation (|2.5I) allows the calculation of any A — > j*N 
decay, once a model for the 7* TV — > A form factors in 
the timelike region is provided. Note however that in 
Eq. (I2.7P the form factors can be directly measured only 
for W = A/a- Consequently, any estimation of the func- 
tion T 1 »N(q;W) has to be done using models that can 
be constrained only in the limit W = Ma- The implica- 
tion is that such models should be largely tested for their 
predictions in other different conditions. Another detail 
on Eq. ([23)1 is that y_ vanishes for q 2 = (W - M) 2 . As 
we discuss later, this point corresponds also to the upper 
limit allowed to q 2 for the reaction to occur. 



B. Explicit expressions for r 7 jv(W) and T e 



T (W) 



We present now the expressions for T~ l ]y(W) and 
r e + e - n(W). The first function is given by Eq. (|2.5I) for 
the q 2 = limit [H, H 113 



T yN (W)=T^ N (0;W). 



(2- 



As for T e + e ~ N (W) it will be determined by integrating 
the function 



r' 



+ e~N 



dT, 



+ e~N 



dq 



(2.9) 



according with 

r e+e - N (W)= T' e+e - N (q,W)dq (2.10) 

J2m e 

Note that the integration holds for the interval 4m;? < 
q 2 < (W — M) 2 , where m e is the electron mass. In this 
case the lowest squared momentum corresponds to q 2 = 
Am 2 , the minimum possible value for a physical e + e~ 
pair. The upper limit is determined by the maximum 
value for q 2 needed for the A with mass W to decay into 
a nucleon (mass M) and is discussed in Appendix \K\ 
The function T' + _ N (q, W) can be determined [IjJ, |3l| 

by 



2a 

T' e+e - N {q,W) = —Y 1 , N {q ] W). 



(2.11) 



The function r' e+e _ N (q, W) diverges when q — > 0, due 
to the presence of 1/q. However, this is not a problem, 



since in (|2.10j) the lower limit of the integration variable 
q (given by 2m e ) prevents the integral from diverging. 

To proceed from here, the calculation of the partial 
widths T lN (W) and r e + e - jv(W) requires a model for the 
7* N — > A form factors in the timelike region, for an arbi- 
trary A mass W. In the past at least three models were 
proposed to this reaction: constant form factors [Til H3 > 
a two-component quark model (model with valence and 
pion cloud components) 0, j34|-l37j and a vector meson 
dominance model from Ref . 39] . In the next section we 
propose a new model based on the spectator formalism. 
This model can be described also as a two-component 
quark model. What is specific of our model is that, in 
addition to the constraints from the spacelike physical 
data, our model was also constrained by the spacelike 
lattice QCD data [3 . 



III. SPECTATOR QUARK MODEL 



We will focus now in the covarian t sp ectator quark 
model for the 7* A -> A reaction [Hf- Here, we will 
describe briefly the properties of the model and summa- 
rize the important results. In its simple version, when 
the nucleon and A are both approximated by an S-state 
configuration for the quark-diquark system, the transi- 
tion form factors are restricted to the dominant magnetic 
dipole form factor, and is decomposed @ into 

G* M (Q 2 ; W) - Gfj(Q 2 ; W) + GljiQ 2 ; W), (3.1) 

where G^ is the contribution of the quark core and GJ# 
represents the effect of the pion cloud. In the previous 
equation W replaces Ma , the physical A mass used in the 
previous applications [(| 0] • Because our original model 
and formulas were developed in the spacelike region, we 
maintain here the use of the variable Q 2 which stands for 
— q 2 . Explicitly Gf^ is written as Q 

G B M {Q 2 ;W) = ^- ] ^L^f v (Q 2 )l{Q 2 ), (3.2) 

where 



1{Q 2 



(3.3) 



is the overlap integral of the nucleon and A radial wave 
functions which depend on the nucleon (P-), the A (P+) 
and intermediate diquark (fc) momenta. The integration 
sign indicates the covariant integration in the diquark 
momentum k; f k = J m^2BB ' wnere Er> = \J m 2 D + k 2 
is the diquark energy (m d is the diquark mass) . Explicit 
expressions for the nucleon radial wave function ip^ and 
the A radial wave function "0a will be presented later. 
As for the factor f v it is represented by 



fv(Q 2 ) = h-(Q 2 ) + 



W + M 
2M 



h-{Q 2 



(3.4) 
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where (i = 1,2) arc the quark (isovector) form 
factors, that parameterize the electromagnetic photon- 
quark coupling [5J, [6|, [8j . The fi± parameterizations will 
be discussed in more detail in the next subsection. 

The pion cloud parametrization GJj was established in 
the physical regime using the factorization Q 



Glf(Q 2 ; W) = 3A 7r G_ D (Q 2 ) 



A 2 



A 2 +Q 2 



(3.5) 



where G D = ^1 + , with Q 2 in GeV 2 , has the 

usual dipole functional form, and A^ and are pa- 
rameters that define the strength and the falloff of the 
pion cloud effects. In particular we take A^ = 0.441 and 
A 2 = 1.53 GeV 2 following Refs. 0, El- More details 
of the model in the physical regime (W = Ma) can be 
found in Refs. 0, 0|- Since the pion cloud parameteri- 
zation given by the right-hand-side of Eq. (|3.5p has no 
explicit dependence on Ma, in its extension to W ^ Ma 
we consider no explicit dependence on W either. Then, 
for W ^ Ma our choice was to keep G* M independent of 
W, that is G^(Q 2 ; W) = G^(Q 2 ; M A ), and the variable 
W could have been dropped in the Eq. (|3.5I) . 

A general comment about the decomposition (|3.ip is in 
order. In the spectator framework the component Gf.j, 



given by Eq. (|3.2[) . is limited by the condition 1(0) < 1, 
which follows from the normalization of the nucleon and 
A radial wave functions and the Cauchy-Schwartz-Holdcr 
inequality. This implies that Gf f (0,M A ) < 2.07 @. 
Since the experimental value is G* M (Q, Ma) — 3 it fol- 
lows that the description of the reaction near Q 2 = is 
not possible, unless the contribution of the pion cloud is 
significant: more than 30% of the total result. The under- 
estimation of G* M (0, Ma) is a result common to several 
models based on constituent quark degrees of freedom 
alone |Q|. 



deep inelastic scattering Q. In the applications we take 
m v — m p (~ m u ), to include the physics associated with 
the p-meson, and Mh = 2M (twice the nucleon mass) for 
effects of meson resonances with a larger mass than the 
p. Note that both functions /j_ and /2_ have a pole at 
q 2 = m p and at q 2 = M 2 . Hereafter we will refer to these 
poles as p-poles and Mh poles, respectively. 

The parametrization (|3. 61) is particularly useful for 
applications of the model to the lattice QCD space- 
like regime. In fact, the decomposition of the current 
into contributions from the vector meson poles (ra p and 
Mh = 2M) is very convenient for a extension of the model 
to a regime where those poles can be replaced by the m p 
and M values given by the lattice calculations, without 
introducing any additional parameters. Examples of suc- 
cessful applications to the lattice regime can be found in 
Refs. [1, El El El- In Refs. [llEji, in particular, one 
can see how well the model describes the lattice data 
from Ref. [II for the j*N — > A reaction, particularly 
for pion masses m n > 400 MeV where the pion cloud 
effects GJjf are suppressed. The valence quark contribu- 
tion 0, EH is also compatible with the estimation of the 
bare contribution from the EBAC model [4lJ. The suc- 
cessful description of the G* M lattice data shows that the 
valence quark calibration of our model is under control. 

To stress the first problem of the extension of (|3.6[) to 
the case q 2 > 0, in Eq. (|3.6[) . we used explicitly the vari- 
able q 2 instead of the variable Q 2 employed in Refs. 0- 
7]: singularities appear at q 2 = m 2 and q 2 = M 2 . The 
larger poles are not problematic for moderated W, since 
as shown in Appendix [2] q 2 < (W — M) 2 . But the 
case q 2 — to 2 has to be taken with care. Such pole is a 
consequence of having the p meson as a stable particle, 
with a zero mass width. One can overcome this limita- 
tion by introducing a finite width T p in the p-propagator 
m 2 /(m 2 — q 2 ), with the replacement m p — > m p 

leads then to the substitution 



non-zero width T 



2 1 p- 



A. Quark current 

In the spectator quark model the electromagnetic in- 
teraction with the quarks is represented in terms of Dirac 
and Pauli electromagnetic form factors, fi± (Dirac) and 
ji± respectively, for the quarks H, [E Hi • Using the vector 
meson dominance (VMD) mechanism, those form factors 
are parametrized as 



f 1± (q 2 ) = X q + (l-X q ) 



2 n 1 



1) 2 



c±- 



M£q 



ml - 9 2 = (M 2 - q 2 ) 



h±{q ) = n±<d± 



m* — q A 



(l-d±) 



M 2 h 



Mf - q 2 



(3.6) 



where m v is a light vector meson mass, Mh is a mass of 
an effective heavy vector meson, k± are quark anomalous 
magnetic moments, c±,d± are mixture coefficients and 
X q a parameter related with the quark density number in 



m- 



m-1. 



ml — q 2 



- q 2 - im p T p 



2 [(m 2 - q 2 ) + impTp] 
(m 2 p — q 2 ) 2 + m 2 p T 2 p 



(3.7) 



Note that this procedure induces an imaginary part in 
the bare quark contributions for the form factors. The 
p-width T p is in fact a real function of q 2 defined only 
for q 2 > 0, as we discuss next, and therefore the results 
in the spacelike regime are unaffected by the redefinition 

The p width can be measured only for the physical 



decay of the p, when q 2 



The experimental value 



is T° = T p (m 2 ) = 0.149 GeV (PDG) [4J. For q 2 > 
one has to consider some parametrization for T p (q 2 ). An 
usual parametrization is |43l - l45l |: 



r P (g 2 



- Ami 



2 \ 3/2 



Am 2 J 



% 2 -4m 2 ), (3.8) 
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where m^ is the pion mass and 9(x) the Heaviside step 
function that cuts the contributions for q 2 < 4m 2 , below 
the 2ir creation threshold (decay p — > 2ir). The previous 
formula includes then the creation of irir states from an 
off-mass-shell p. Equation (|3 . 8|) assures that there is no 
width near q 2 = 0. Therefore the imaginary contribution 
appears only for q 2 > 4to 2 ~ 0.076 GeV 2 . 



B. Scalar wave functions 

The radial (or scalar) wave functions taken in this 
work, respectively for the nucleon and A, are 

N N 



TpN(P,k) 



m D (/3i + Xat)(/?2 + Xn) 
N A 



(3.9) 
(3.10) 



m D (a + xa) 3 ' 

where mp is the diquark mass, /?2 and a are momen- 
tum range parameters (in units mp) and 

(Mb - mp) 2 - (P - k) 2 

Xb = 77 , (3.11) 

Mb tod 

for B = N (M B = M) and B = A (M B = Ma), is a vari- 
able without dimensions that includes the dependence in 
the quark momentum (P — k) 2 . As for Nb, (B = N,A 



they are positive normalization constants. See Refs. j^. 
for details. The representation of the wave function in 
terms of xb given by Eq. (|3.1ip has advantages in the 
applications to the lattice regime @, [IH, [l]| . 

The scalar wave functions are important for the present 
calculations because they are part of the overlap integral 
defined by Eq. p.3[) . To apply the expressions to the 
timclikc region one has to choose a configuration with 
Q 2 < (q 2 > 0). That can be achieved by considering the 
reaction 7* A — > A in the A rest frame, with the following 
configuration: P + = (W, 0) as the A momentum and 
P_ = (En, — q), with En = \J M 2 + q 2 , as the nucleon 
momentum. In those conditions the photon momentum 
is represented by q = P + — P-, as q — (w, q), where 

W 2 - M 2 + q 2 



2W 
(W 2 + M 2 



q 



2 \2 



AW 2 



M 



(3.12) 



Those variables correspond to the timelike region when 
< q 2 < (W - M) 2 . See details in Appendix EJ 

Using Eq. (|3.2I) and the integral (|3.3I) for the kine- 
matics (|3.12j) . together with the extension of the current 
given by (|3.7[) , one can calculate the contribution for Gfj . 
Note that as the function (|3.7I) has an imaginary compo- 
nent, G?, is now complex. 



C. Pion cloud contribution 

The most phenomenological part of the model pre- 
sented here is the parametrization of the pion cloud con- 



tribution through Eq. (|3. 51) . Although the valence quark 
parametrization has been validated by lattice QCD sim- 
ulations and the EBAC estimations of the quark core 
contributions [13, EH, the contributions from the pion 
cloud were estimated only phenomenologically. In fact 
they were extracted directly from the physical data, af- 
ter the calibration of the valence quark effects [H, 0] • 

For the pion cloud component of the form factor we will 
compare two different generalizations of Eq. (|3.5[) for the 
timelike region. We start with a simple model, a naive 
generalization of the model from Refs. H El t0 thc 
timclike region. Next we discuss the possible limitations 
of that approach and introduce a different parametriza- 
tion motivated by the expressions for the pion cloud de- 
rived from xPT. 



1. Naive model (model 1) 

In a first approach we took the pion cloud contributions 
for the G* M form factor by Eq. (|3.5[) . as in the spacelike 
regime, but now evaluated in the timelike kinematic re- 
gion. We have to take into consideration now the poles 
for q 2 > (Q 2 < 0). We re-write Gd as 



G D (q 2 ) = 



A 2 

D 



A 2 



T 



(3.13) 



where K 2 D 



0.71 GeV 2 is the cutoff of the dipole form 
factor. As it happens to the p— term in the quark current, 
also this factor has a pole at q 2 = 0.71 GeV 2 , but in this 
case it is a double pole. We apply the procedure used 
before to the p propagator, i.e. definiting a width to the 
function Gd, by making 



G D (q 2 



A 2 



(A 2 D -q 2 ) 2 



A 2 r 2 

D D 



(3.14) 



[(A 2 D -q 2 ) 2 -A 2 D T 2 D+l 2(A 2 D -q 2 )A D r D ] 



where r_o is thc width associated with its pole 



poles m 2 and A 2 D are close (q 



0.6 GeV 2 versus q 2 



As the 
0.7 

GeV 2 ), we will use T D (q 2 ) = T p (q 2 ). With G D defined 
as above, also G\ { is a complex function in the timelike 
regime. 

As for the extra dipole factor in Eq. Q3.5P : ( A % 

where in the applications A 2 ~ 1.5 GeV 2 , it is far way 
from the p— poles region. For A masses not very large 
compared with Ma the possible effect of the finite width 
is less significant since q 2 < (W — M) 2 < A 2 . 

We call the model defined by Eqs. ((331) and (|3"J4]) 
model 1. The result of the extension of the model to 
the timelike region for the case W — M A and T p = 
is presented in Fig. Q] We used the parametrization of 
Ref. [1 2| | for the valence quark contributions, but ne- 
glected the D-state contributions (< 1%). As in this 
case Q 2 > —(Ma - M) 2 ~ -0.086 GeV 2 , and therefore 
the allowed values for Q 2 are far way from the Q 2 < 



fi 




0.4 0.6 

Q 2 (GeV 2 ) 



FIG. 1: 7* TV — Y A(1232) magnetic form factor in the timclikc 
and spacelike region given by the model of Ref. [T3| . The data 
for Q 2 > is the same as the one presented in Ref. [12| . The 
data for Q 2 = is from Ref. [46| , The dashed line represents the 
contributions of the valence quark core (bare contribution). The 
solid line is the result of the sum of valence quarks and pion cloud 
contributions. 



poles, the corrections due to the imaginary components 
are small. In the figure we show also physical data for 
Q 2 > and the result for Q 2 = from Ref. 0. 



2. xPT motivated model (model 2) 



Instead of Eq. (|3.5p for the pion cloud effect we can 
use a different parametrization, based in a different com- 
bination of multipolc functions. For instance, in the two- 
component model from Refs. [l9|, [37| the contribution 
from the pion cloud is proportional to the function F p , 
interpreted as the p— propagator, derived from xPT. This 
function F p was presented in Refs. [351. 1371 l38j taking into 
account the pion loop contributions to the p— propagator. 
Here we simplified the exact expression in those refer- 
ences by assuming its limit when q 2 3> 4mJ, and using 
the normalization F p (0) = 1. For Q 2 = —q 2 > we 
obtained then 



F P (q 2 ) 



(3.15) 



In the previous equation, the physical p-width F° was 
taken to be T° p = 0.149 GeV. Reference pjj], uses instead 
T° p = 0.112 GeV. 

Equation ()3.15j) . derived in the low q 2 chiral pertur- 
bation regime, has a faster falloff for the p-propagator 
[with l/(Q 2 log^-)] than model 1 [with 1/Q 2 ] for large 

Q 2 . We used it here to explore alternative parametriza- 
tions to model 1 for the pion cloud contributions. The 
parametrization for the pion cloud contribution (|3.5[) is 

/ A 2 

proportional to the dipole factor I ^ * c 




0.2 



0.4 0.6 
Q 2 (GeV 2 ) 



0.8 



FIG. 2: Comparing the dipole form factor with the function 
Fp{Q 2 ) given by Eq. pT5]l , 



is a large cutoff, and also to Gn, the dipole form fac- 
tor. Although the dipole factor depending on A„. was 
chosen phenomenologically and determined by a fit to 
the data, one has no reason a priori to use the particu- 
lar form of Gd to parameterize an extra falloff 1 of G\ L . 
The inclusion of Gd was motivated by the traditional 
convention of dividing the form factor G* M by 3Gd when 
showing results. With the parametrization (|3.5|) one has 
an asymptotic dependence of G\. s oc 1/Q 8 . 

We note that in the spacelike region where the pion 
cloud effects are more important, < Q 2 < 1 GeV 2 , the 
functions Gd and F p give very similar results, as seen in 
Fig. [21 This suggests that one can also use 



GUQ^^iKF^q 2 ) 



A 



Al-q 2 



with 



F p {q 2 ) 



2 _ lH q 2 log JJ^ +l H q 2' 



m 2 - q 



(3.16) 



(3.17) 



to extend Eq. (|3.15l) to the timelike kinematics. 2 The 
imaginary part in Eq. (|3.17j) is a consequence of two pion 
production (or transition p 2ir) which is possible in 
the timelike region when q 2 > 4m 2 . 

We will call the model defined by Eq. (|3.16p model 
2. One implication of the new form for the pion cloud 



where A n 



The function G r> provides a good approximation for the behav- 
ior of the nucleon electromagnetic form factor at low Q 2 . 
2 In the transformation from Eq. I|3.15p to Eq. I|3.17p there is an 
ambiguity from the factor 

log(-l) = i-K + i(2n)t, 

where t is an integer. In this case the ambiguity is fixed by the 
sign of the imaginary part from Ref. Il9|| . 
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FIG. 3: Results for the modulus of the ■yN — > A mag netic form factor for different A masses (W). At left: model 1; at right: model 2. 
In both cases the valence quark model is from Ref. [12 ] . 



contributions is a falloff as l/(q 6 logq 2 ), slower than for 
model 1 [l/q s falloff]. [Note that G D oc l/q 4 and F p cx 
l/(<Z 2 logg 2 )] 

Another important feature of the function (|3.16[) is 
that its imaginary part does not peak for q 2 = m 2 ~ 0.6 
GeV 2 but for q 2 ~ 0.3 GeV 2 , because of the logarithm 
corrections. That effect changes the Q 2 dependence of 
the pion cloud contributions to G* M , relatively to model 
1. 

A note about F p (q 2 ) given by Eq ( 3.171): i t was derived 
from the exact result in Refs. [TiJ 35, 38] in the limit 
q 2 ^> 4m 2 ~ 0.08 GeV 2 , as explained before. However, 
we have checked that our simplified formula, although 
approximate, does not deviate too much from the exact 
one, even when that limit does not hold. Therefore our 
formula gives a good qualitative description of the chiral 
behavior in the whole domain q 2 > 0. We note only that 
the results from Eq. (|3.17p and the results in Refs. [l9|, 
l35l [381 ] differ in a slight deviation of the location of the 
peak of the imaginary part of F p . Using Eq. (13.171) the 
peak is at q 2 ~ 0.3 GeV 2 , while in Ref. [H, HEHg the 
peak is at q 2 ~ 0.4 GeV 2 . 

Finally, as for the quark current (|3.6[) in the bare quark 
contributions, we will not replace our parametrization of 
the p-propagator (|3.7p by (|3.15p since they differ substan- 
tially and our parametrization was alreadycalibrated by 
the physical data [1-0] and lattice data [H, [H[ for the 
nucleon and A(1232) systems, in the spacelike region. A 
different parametrization for timelike and spacelike re- 
gions would be inconsistent. 



IV. RESULTS 

We will divide the presentation of our results into two 
parts. In the first part we show the results of the mag- 
netic form factor G* M , calculated with the two models 
described in the previous section. In the second part we 



show the results for the width functions T 7 n(W) and 
T e + e ~ N (W), and also for the A mass distribution func- 
tion g&(W). 



A. Form factors in the timelike region 

Contrarily to what happens in the spacelike domain, 
in the timelike region the form factor G* M has a non-zero 
imaginary part. Because of Eq. (12.71) we are interested in 
the absolute value of the form factor, \G* M |, which enters 
into |Gt(<? 2 ; W)\ 2 . Although the form factors are defined 
for any value of W > M, we show here results for selected 
values of W only. We recall that the range of Q 2 for the 
function G* M (Q 2 ;W) depends on W, as established by 
the condition Q 2 > ~(W - M) 2 . 

The results for \G* M \ at energies W = 1.232, 1.500 and 
1.800 GeV, for both models, are presented in Fig.[3l One 
notes that the value of \G* M \ near Q 2 = decreases with 
W. We will see that this is a consequence of the valence 
quark contribution given by Eq. (|3.2p . The same effect 
was observed in lattice QCD simulations where large pion 
masses induce large nucleon and A masses [l2|,|4(]|]. In the 
figure it is also clear that the two models differ substan- 
tially in the Q 2 dependence of |Gj^|. For model 1, |G^| 
increases as Q 2 decreases, for all values of W, and this 
behavior is enhanced as W becomes larger. A peak (not 
shown in the graph because it is too large) is present, 
near q 2 = A 2 D ~ 0.71 GeV 2 when W = 1.800 GeV. A 
first conclusion is therefore that model 1 generates very 
strong, and probably unphysical contributions to |G^| 
in the timelike region. Model 2, in contrast, gives mod- 
erated contributions only (larger than in the spacelike 
region but with the same magnitude) and is therefore a 
much more reasonable model. To better compare the two 
models we have to analyze the real and imaginary parts 
of G* M separately. 

Since the valence quark contributions are common to 
both models, we start by looking to the bare term Gf f . 
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FIG. 4: Valence contribution for -yN — > A(1232) magnetic form 
factor for W = 1.232. 1.500 and 1.800 GeV. 



The results are presented in Fig. 2) For Q 2 « and val- 
ues of W not too large when compared with W = 1.232 
GeV, the real part dominates, as expected from the re- 
sults for the physical case (W = Ma). For larger values 
of W and low Q 2 the real part dominates increasingly 
less. As for the imaginary part, we recall that is zero 
down to Q 2 ~ -0.08 GeV 2 (because T p = 0). But as 
Q 2 decreases, for larger W values, we can observe the 
effect of the p-mass poles emerging at Q 2 = 



2 ~ -0.6 

GeV 2 , and strengthening the imaginary parts of G^ . For 
W = 1.800 GeV the real part also increases in the re- 
gion Q 2 < —m 2 due to the impact of the p-mass poles. 
In this case, however, the other terms from the VMD 



parametrization (|3.6[) , the constant term and the Mh- 
mass poles, are relevant as well. All these contributions 
to G£f are balanced, and therefore reduced, in the fi- 
nal result, by the pion contributions G* M to be discussed 
next. 

We turn now to the term G\ f which gives the pion 
cloud contribution, and where the two models differ in 
the timclikc region. The results are presented in Fig. [5] 
for the same values of W as before. Comparing the two 
models, we can say that both have similar results for the 
real part of G\[ in the region —0.2 GeV 2 < Q 2 , but differ 
significantly for smaller values of Q 2 (larger values of q 2 ). 
That is the consequence of the double pole q 2 = K 2 D in 
the pion cloud formula for model 1. In model 2 there is 
no such contribution from the pion cloud and the values 
for real and imaginary parts are more moderate. Note 
that the strong peak for W = 1.800 GeV at Q 2 ~ -0.65 
GeV 2 , for model 1 is a consequence of Eq. (|3.14p . and 
differs from model 2 by an order of magnitude. 

We look now to the imaginary part of Gf.j. Our first 
observation goes to the imaginary part of G\ t in the re- 
gion near Q 2 = 0. In model 1 it is identically zero for 
-0.076 < Q 2 < [because T p = from Eq. but in 

model 2 it is different from zero, although small, (this is 
a consequence of the approximation considered in func- 
tion F p discussed previously). The second observation 
is that the significant difference between models 1 and 2 
is the sign of the imaginary part of G\f. model 1 gives 
positive contributions, while model 2 gives negative con- 
tributions. This model is motivated by xPT, satisfies 
chiral constraints for the p propagator, which has a non- 
analytical pole near Q 2 ~ —0.3 GeV 2 present in F p and 
with its origin in the pion loop contributions. 

With this detailed analysis we come to understand the 
large difference between the two models shown in Fig. [3] 
The figure provides a strong indication that model 1 is 
not a reasonable model: the results from this model are 
strongly dominated by the pole q 2 — K 2 D — 0.71 GeV 2 
which induce extremely large (and probably unphysical) 
contributions for increasingly large W values. On the 
other hand, model 2 contains the input from x?T, and 
has therefore a more solid basis. The results are very 
sensitive to this input and clearly exclude model 1. 

From the previous discussion we favor the results from 
model 2. The final results (bare quark core plus pion 
cloud) for both the real and imaginary part of the form 
factor G* M are presented in Fig. [5] One realizes that 
the real part dominates over the imaginary part for 
Q 2 > —0.15 GeV 2 . We can say that the dominant con- 
tributions to the imaginary part are the poles from GJ^ 
(induced by F p ) around Q 2 = —0.3 GeV 2 , and the p- 
terms in the quark current from the bare contribution. 
The effect of those poles is particularly evident in the 
curve for W = 1.800 GeV. 

The main and general conclusion from Figs. 21 [5] and 
[6] is that the final structure for G* M emerges from a com- 
bination of effects, namely from the VMD model poles 
and also the pion cloud effects. The two processes in- 
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FIG. 5: Pion cloud contributions (G^ f ) for the G* M form factor (real and imaginary parts). The results from models 1 and 2 are shown. 



terfere crucially and determine the structure of the hnal 
amplitude. 

Later, we will discuss the applicability of the model 
for W > 1.8 GeV, the effect of the remaining poles, and 
the impact in the observables in consideration. We em- 
phasize that any extension of the A form factors to the 
timelike region has to rely on models and cannot be di- 
rectly estimated from experimental data only. It is then 
important to compare our model with models with a sim- 
ilar content, such as the two-component quark model of 
Ref. [l9j, also defined in the timelike region. In this last 
model the contribution from the coupling to the quark 
core (valence contribution) is 0.3% near Q 2 = (99.7% 



of pion cloud), while in our model one has 55.9% (44.1% 
of pion cloud). This significant difference between the 
contributions of the quark core is due to a different, and 
somewhat arbitrary, classification of the two effects. In 
the model of Ref. [19l [35| the term from the pion cloud 
is also classified as an effective part of the VMD mecha- 
nism, since it is proportional to the function F p for the p 
propagator. Therefore, in that model the VMD mecha- 
nism/pion cloud term is the only relevant effect fl9l l35l]. 
In our formalism, the coupling with the quarks is cali- 
brated directly by a VMD parametrization and although 
it gives the dominant contribution, it is not the only one 
to affect the results. Our model has the advantage of 
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FIG. 6: Real and imaginary components of the form factor G* M 
for model 2. 



having been tested successfully by the lattice QCD simu- 
lations (in a regime where the pion cloud is small), and of 
agreeing with the EBAC data analysis for the bare quark 
core contributions to the pion photoproduction data [l2| . 
These tests suggest that our estimation of the quark core 
structure is under control, since the model is largely con- 
strained in a variety of kinematic domains. Another im- 
portant point is that our model allows a direct physical 
interpretation of the parameters involved, in terms of the 
range of the baryon wave functions. 



B. Results for Y lN {W) and T e+e - N (W) 

We will discuss now the partial widths T 7 n(W) and 
r e + e - n{W). We will also show T n ^(W), given by 
Eq. (12.31) . together with the calculation of <7a {W), de- 
fined by Eq. ([2~T]) . 

We start by showing in Fig. [7] the function 
dr ^f N {W;q) for the cases W = 1.232, 1.500 and 1.800 
GeV. This figure includes the results from model 1 
(dashed line), model 2 (solid line) and also the result 
of a calculation where the form factor G* M is taken as 
constant, defined by the value of G* M (W, 0) at the pole 
W = Ma (dotted line), given by the experimental value 



G* M (Q 2 



G 



3.0]. This last case was also 
considered in Ref. [19( and it is useful as a reference for 
the q 2 dependence of our results. The figure illustrates 
that, in line with the results in the previous subsection, 
for model 1 r' e+e _ N (q, W) is enhanced for large q and 
large W values (see result for W — 1.800 GeV). 

To determine the di-lepton production width, 
r e + e - jv(W), one has to integrate Eq. (|2.11l) using 
Eq. (I2.10[) . This is equivalent to calculate the integral of 
the functions represented in Fig. [7] for each value of W 
in the interval [2m e , W-M]. Therefore T e+e - N (W) = 
when W < M + 2m e . The calculation of the function 
T 7N (W) proceeds through Eq. (12~B . The results 
obtained for the two widths within the three models 
discussed before, are in Fig. [5] 

Finally, T^nIW) is estimated using Eq. (|2.3j) and the 
function 



y/[(W + MY - ml] [(W - M) 2 - 
2W 



',(4.1) 



defined for W > M + m v and q-n{W) = otherwise. 
^ttn(W) is then a positive function for W > M + m n . In 
Fig. [9] we present the three partial widths obtained with 
model 2, the one that we favor for the reasons explained 
in the previous subsection. 

We turn now to the A mass distribution function 
5a(W) defined by Eq. (|2~Tj) . As the channel A -> 7^ 
is largely dominant, T to t{W) ~ Ykn{W) and the nor- 
malization of <7a(W0 can be done in that approxima- 
tion. Considering T^ N (M A ) ~ r tot (M A ) ~ T e ^ p , with 
the experimental result Y e £ p ~ 0.118 GeV [42j], one has 
A = 0.7199. 

The results for the partial contributions to gA (W) are 
given by 



T^ N {W) 

Ttot{W) 



9a(W) 
n(W) 



T tot (W) 



g A (W), (4.2) 



and are shown in Fig. [TU] for the constant form factor 
model (dotted line), model 1 (dashed line) and model 2 
(solid line). The total results for <?a(W) are shown in 
Fig. El for the model 2. 
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FIG. 8: Partial widths as function of W for different models. 
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Note that models 1 and 2 have the same result for the real decay 
A — > jN, because the models are identical at Q 2 = 0. 



We restricted our results to the region W < 2 GeV. 
Above that region one has to take into account the addi- 
tional pole structure of the form factor G* M that appears 
for large q 2 values. Another reason for not having our 
model applied to larger W values is that, for W > 1.5 
GeV, the reactions are expected to be dominated by res- 
onances as the TV* (1440), AT* (1535), A(1600) among oth- 
ers, instead of the A(1232) alone. 

As for the first problem, one has to find a way deal 
with those large q 2 singularities. From Eq. (|3.16p for 
the pion cloud component, already, the pole at q 2 — 
A 2 ~ 1.53 GeV 2 , makes the form factor diverge for 
W > M + A w ~ 2.17 GeV, since q 2 < (W - Mf . Also, 
from the valence quark component one has a singular- 
ity q 2 = AT 2 ~ 3.53 GeV 2 . To remove those singulari- 
ties we may introduce an effective width Tx, in analogy 
with Eq. (|3.7p or (|3 . 14[) . for single or double poles, with 
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FIG. 9: Partial widths as function of W for model 2. The total 
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FIG. 10: Partial contributions of the channels A — > 7./V and 
A — > e + e~ N (in GeV -1 units). The constant form factor model is 
represented by the doted line; the model 1 by the dashed line and 
the model 2 by the solid line. Note that models 1 and 2 have the 
same result for the real decay A — > 7JV, because the models are 
identical at Q 2 = 0. 



a constant width T x . This procedure adds a new pa- 
rameter to the model. We started by verifying that a 
large width T x (for instance T x ~ 1 GeV) will affect the 
results for W < 2 GeV, while a very small width will in- 
duce significant oscillations in the functions gA^yN(W) 
and #A-s.e+e- n(W) (enhancement of F' e + e - N (q, W) near 
the poles). By taking T x — 2r° the results obtained 
for W < 2 GeV are almost unchanged, and the results 
for W > 2 GeV are also smooth functions of W. With 
this choice the smooth dependence on W of the functions 
<7a^7At(W) and 3A^e+e-Af(^)i obtained with model 2 
is maintained for higher W values. We stress, neverthe- 
less, that the application of our model to large W values 
is to be viewed only as an extreme test to its limits. 



V. CONCLUSIONS 

In this work we have presented a covariant approach 
to describe the the Dalitz decays A — > jN, A — > 7* A, 
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and the sum off all decay channels (A — > nN, A — > ■yN and A — > 
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and we have calculated the A mass distribution function 
g&(W) within that approach. Our framework can be 
used to simulate the NN — > e + e~ NN reaction at mod- 
erate beam kinetic energies (1 — 2 GeV). The code used 
in this work can be supplied under request. 3 

Our calculations are based on an unified description of 
the 7* A — > A reaction, in both the spacelike and time- 
like regimes. We start with a model tested previously 
jacelike physical and lattice QCD simulation data 
and generalize it to the timelike regime. In this 
formalism the electromagnetic interaction can be decom- 
posed into two mechanisms: the direct photon coupling 
with the quarks and the interaction with the pion cloud. 

For the first mechanism we extended the quark cur- 
rent and wave functions, obtained for q 2 < 0, to the 
region q 2 > without additional changes, except for a 
non-zero width of the effective vector mesons included in 
the VMD parametrization of our quark electromagnetic 
current. For the pion cloud contribution we probed two 
different parametrizations: a naive generalization of our 
spacelike model and a more elaborated model based on 
xPT. Although the two models behave very similarly in 
the spacelike region, they differ substantially in the time- 
like region. 

The results of the models for the form factor G* M , as 
are used to calculate the par- 
n{W) and the partial mass 
n{W). 



a function of q and W 
tial widths T lN (W), r A ^ e+ 
distributions functions gA^~/N{W), g& 

A first important conclusion of this work is that the q 2 
dependence of the form factor has an impact in the final 
results, and has therefore to be under control: the results 
for the A mass distribution functions, where the form 
factor is taken with its full q 2 dependence, can differ by 
a factor of 4 from the results obtained with the constant 



form factor model. 

We verified also that the results are sensitive to the 
analytical extension of the pion cloud parametrization 
to the timelike region. Model 1 (naive model) generates 
unreasonably large contributions for moderate squared 
momentum q 2 , for larger values of W, as consequence of 
a spurious (not physically motivated) pole q 2 = K 2 D in the 
pion cloud contribution. Model 2 is motivated by X-PT 
and includes the function F p for the p propagator with a 
non-algebraic pole near q 2 ~ 0.3 GeV 2 originated by pion 
loop corrections. This model gives smooth contributions 
to the partial A mass distribution functions gA^yN{W) 
and 5A^-e+ e -Af(W / ) that vary slowly with W for large W. 

A second important conclusion is then that our calcu- 
lations support the need to have under control the effect 
of the pion cloud contributions. Our framework is suit- 
able for this because its bare quark core component is 
constrained by experimental data and lattice QCD sim- 
ulation data. In addition, the fact that the pion cloud 
content of model 2 is consistent with xPT makes model 2 
reliable, at least in its domain of validity. The results of 
model 2 give moderated contributions for |G^f| 2 (since 
it has no spurious pole at q 2 — A|>), which are deter- 
mined by the combined effect of two important features: 
the pion cloud term structure near q 2 = 0.3 GeV 2 , and 
the quark core contributions from the p— pole included in 
the VMD structure of the quark current, near q 2 = 0.6 
GeV 2 . 

We may say that while spacelike data does not con- 
strain models sufficiently well enough, timelike data for 
dr "dq N f° r different values of W (see Fig. [7]) are impor- 
tant and necessary to select between models in a decisive 
way. In the A case discussed in this work this is spe- 
cially true for the pion cloud effects. But the timelike 
data can also be useful to calibrate the widths and high 
mass poles of the VMD parametrization of the current 
needed in valence quark component, particularly for res- 
onances heavier that the A(1232). 

For high values of W the assumption that the A(1232) 
resonance is the only state playing a role in the reac- 
tions becomes questionable. In the regime W > 1.5 GeV 
other resonances can be relevant, as the spin 1/2 reso- 
nances A*(1440) and A*(1535). Once those states are 
calibrated for the Q 2 = - 2 GeV 2 region, one can ex- 
tend the models from Refs. [HI, [l6| to the timelike region 
too. It is also expected that the spin 3/2 channels as the 
A(1600) state are important for large W. This defines a 
study of high interest, since the constraints in the space- 
like region are very scarce, and the available data at the 
photon point suggest a strong contribution from the pion 
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Future applications of our formalism can include the 
study of the 7* A — > N reaction where the final nucleon 
has an arbitrary mass W, also in the timelike region. 
Since the quark current was already defined in the time- 
like region and we have already a model for the nucleon 
system || , no additional ingredients are necessary. Such 
study may provide an important theoretical input to the 
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study of the reaction pp 



e in complement to the 
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then q as the photon three momentum (symmetric to the 
nucleon three momentum) in the final state, and w as the 
photon energy in the R rest frame. All those variables 
are related with W 2 and q 2 according to 



2\2 



(W 2 +M 2 - q 2 ) 
AW 2 
Pr ■ q W 2 



M 2 



M 

„2 



w 



2W 



2W 



(A2) 
(A3) 

From 



In the last relation we used En 
Eq. (TA21) 

creases. As q 2 > 0, it can be proved that there is an 
upper limit to q 2 (given by the condition q 2 =0). 
upper limit is then 



we can conclude that q 2 decreases when q 2 in- 



The 



Appendix A: Kinematics 

We consider here the final state of the decay process 
of a resonance R according to R — > 7* A. Assuming W 
as the invariant mass of the resonance, we can write the 
four-momentum in the rest frame of R as Pr — (W,0). 
In this frame we can also write 



P N = (E N ,-q), 



q = (w,q), 



(Al) 



where q is the photon momentum: Pr = Pn + q, with 
En = \J M 2 + q 2 (M is the nucleon mass). We define 



q 2 \ =(W-M) 2 . 



(A4) 



This is then the largest value of q 2 for which the timelike 
form factors are defined. 

In conclusion, the timelike form factors are defined only 
in a limited interval [0, (W—M) 2 ]. One has then different 
ranges according to the value of W. In the case W = M 
we have only q 2 = 0. At the pole W = 1.232 GeV the 
maximum value of q 2 is 0.0859 GeV 2 . For W = 1.500 
GeV and W — 1.800 GeV the maximum value of q 2 is 
respectively 0.315 GeV 2 and 0.741 GeV 2 . 
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